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Abstract 



A search for supersymmetry in final states with jets and missing transverse energy is 
performed in pp colHsions at a centre-of-mass energy of = 7 TeV. The data sample 
corresponds to an integrated luminosity of 4.98 fb^^ collected by the CMS experiment 
at the LHC. In this search, a dimensionless kinematic variable, ocj, is used as the main 
discriminator between events with genuine and misreconstructed missing transverse 
energy. The search is performed in a signal region that is binned in the scalar sum 
of the transverse energy of jets and the number of jets identified as originating from 
a bottom quark. No excess of events over the standard model expectation is found. 
Exclusion limits are set in the parameter space of the constrained minimal supersym- 
metric extension of the standard model, and also in simplified models, with a special 
emphasis on compressed spectra and third-generation scenarios. 
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1 Introduction 

Supersymmetry (SUSY) is generally regarded as one of the likely extensions to the standard 
model (SM) of particle physics HHHl. It is based on the unique extension of the space-time sym- 
metry group underpinning the SM, introducing a relationship between fermions and bosons. 
A low-energy realisation of SUSY, e.g. at the TeV scale, is motivated by the cancellation of the 
quadra tically divergent loop corrections to the Higgs boson mass in the SM IZHSl. These correc- 
tions are proportional to the masses of the particles that couple to the Higgs boson. The most 
relevant terms come from the interplay between the masses of the third generation (top and 
bottom) squarks, and the largest Yukawa coupling (of the top quark). 

In order to avoid large cancellations in these loop corrections, the difference in masses between 
the top quark and the third generation squarks must not be too large [9J. While the majority of 
SUSY particles might not be accessible at the present energy and luminosity delivered by the 
Large Hadron Collider (LHC), the recent discovery of a low-mass Higgs boson candidate llTOl 
[Tl] motivates models in which top and bottom squarks appear at the TeV scale. Furthermore, if 
the multiplicative quantum number R-parity [12J is conserved, SUSY particles will be produced 
in pairs and decay to SM particles and the lightest SUSY particle (LSP), which is generally 
assumed to be weakly interacting and massive. This would result in a final state that is rich 
in jets, especially those originating from bottom quarks, and contains a significant amount of 
missing transverse energy, Ut- 

This paper summarises a search that is designed to be sensitive to missing transverse energy 
signatures in events with two or more energetic jets that are categorised according to the num- 
ber of reconstructed jets originating from bottom quarks (b-quark jets) per event. With respect 
to previous searches ||T3l[T4| , this refinement provides improved sensitivity to third generation 
squark signatures. However, the same inclusive search strategy is deployed, thus maintaining 
the ability to identify a wide variety of SUSY event topologies arising from the pair production 
and decay of massive coloured sparticles. 

The ATLAS and Compact Muon Solenoid (CMS) experiments have performed various searches 
IIT3ti211 for the production of massive coloured sparticles and their subsequent decay to a final 
state of jets and missing transverse energy. These searches were performed with a dataset 
of pp collisions at \/s = 7 TeV, and no significant deviations from SM expectations were ob- 
served. The majority of these searches have been interpreted in the context of a specific model 
of SUSY breaking, the constrained minimal supersymmetric extension of the standard model 
(CMSSM) [22- 24 1 . The simplifying assumption of this model is universality at an energy scale 
of 0(10^^) GeV which makes the CMSSM a useful framework to study SUSY phenomenology 
at colliders, and to serve as a benchmark for the performance of experimental searches. 

However, the universality conditions of the CMSSM result in significant restrictions on the 
possible SUSY particle mass spectra and thus kinematic signatures. This limits the interpreta- 
tion of the results in scenarios such as the direct production of third-generation squarks and 
compressed spectra, where the mass difference between the primary produced sparticle (e.g., 
a squark or a gluino) and the LSP is small. Therefore, in order to complement the interpreta- 
tion within the CMSSM, simplified models p5]-l27l are also used to interpret the search results. 
These models are characterised using a limited set of SUSY particles (production and decay) 
and enable comprehensive studies of individual SUSY event topologies. The simplified model 
studies can be performed without limitations on fundamental properties such as decay modes, 
production cross sections, and sparticle masses. A special emphasis is placed on interpretation 
within models involving compressed spectra or third generation squarks. 
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3 Object definitions and event reconstruction 



2 The CMS apparatus 

The central feature of the CMS detector is a superconducting solenoid, which provides an axial 
magnetic field of 3.8 T. The bore of the solenoid is instrumented with several particle detection 
systems. Silicon pixel and strip tracking systems measure charged particle trajectories with full 
azimuthal (^) coverage and a pseudorapidity acceptance of I // 1 < 2.5, where?/ = — ln[tan(0/2)] 
and 6 is the polar angle with respect to the counterclockwise beam direction. The resolutions on 
the transverse momentum (pj) arid impact parameter of a charged particle with pj < 40 GeV 
are typically 1% and 15 f/m, respectively. A lead-tungstate crystal electromagnetic calorimeter 
(ECAL) and a brass /scintillator hadron calorimeter surround the tracking volume. The region 
outside the solenoid is covered by an iron/ quartz-fiber hadron calorimeter. The ECAL covers 
< 3.0 and provides an energy resolution of better than 0.5% for unconverted photons with 
transverse energies above 100 GeV. The hadron calorimeters cover < 5.0 with a resolution 
in jet energy, E (GeV), of about 100%/ \/E for the region |f/ 1 < 3.0. Muons are identified in gas- 
ionization detectors, covering \rj\ < 2.4, embedded in the steel return yoke. The CMS detector 
is nearly hermetic, which allows momentum-balance measurements in the plane transverse to 
the beam axis. A two-tier trigger system is designed to select the most interesting pp collision 
events for use in physics analysis. A detailed description of the CMS detector can be found 
elsewhere p8|. 

3 Object definitions and event reconstruction 

The event reconstruction and selection criteria follow the procedure described in Refs. ||T3l[T4| . 
Jets are reconstructed from energy deposits in the calorimeters, clustered by the anti-Zcj algo- 
rithm p9| with a distance parameter of 0.5. The raw jet energies measured by the calorimeter 
systems are corrected to establish a uniform relative response in tj and a calibrated absolute 
response in transverse momentum with an associated uncertainty between 2% and 4%, de- 
pending on the jet )] and pj II30I . Jets considered in the analysis are required to have transverse 
energy Ej > 50 GeV and the two highest-Ej jets must each satisfy Ej > 100 GeV. These two 
Ej requirements change under special circumstances described in Section |4j The highest- Ej 
jet is additionally required to be within the central tracker acceptance {\rj\ < 2.5). Events are 
vetoed if any additional jet satisfies both Ej > 50 GeV and \r]\ > 3, or rare, spurious signals are 
identified in the calorimeters IST). To suppress SM processes with genuine Kt from neutrinos 
in the final state, an event is vetoed if it contains an isolated electron [32J or muon [33] with 
Pt > 10 GeV. Further, events with an isolated photon [34J with pj > 25 GeV are also vetoed. 

The presence of a b-quark jet is identified through a vertex that is displaced with respect to the 
primary interaction, using an algorithm that attempts to reconstruct a secondary vertex using 
tracks from charged particles associated to each jet. Using a likelihood ratio technique, the 
combined secondary vertex algorithm |35| incorporates several variables related to the vertex, 
such as decay length significance, mass, and track multiplicity, to build a discriminator that 
distinguishes between jets originating from bottom quarks and those from other sources. These 
include jets from charm quarks (c-quark jets) and light-flavour quarks. The algorithm also 
provides a value for this discriminator based on single-track properties, when no secondary 
vertices have been reconstructed. Discriminator values above a certain threshold are used to 
tag jets as originating from b quarks. This threshold is chosen such that the mistagging rate, the 
probability to tag a jet originating from a light-flavour quark, is approximately 1% for jets with 
transverse momenta of 80 GeV P35ll36 l. The same threshold results in a b-tagging efficiency, the 
probability to correctly tag a jet as originating from a bottom quark, in the range 60-70% [i35l 
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The following two variables characterise the visible energy and missing momentum in the 

transverse plane: the scalar sum of the transverse energy Ej of jets, defined as Hj = Ylfj^l £t''/ 
and the magnitude of the vector sum of the transverse momenta pj of jets, defined as Ht = 

I Tjf^l pV' \, where Njet is the number of jets above the Ej threshold. Significant hadronic activity 
in the event is ensured by requiring Hj > 275 GeV. Following these selections, the background 
from multijet production, a manifestation of quantum chromodynamics (QCD), is still several 
orders of magnitude larger than the typical yields expected from a SUSY signal. 

4 Selecting events with missing transverse energy 

The IX.J kinematic variable ||T3ll37| is used to efficiently reject multijet events without significant 
£t/ including those with transverse energy mismeasurements, while retaining a large sensi- 
tivity to new physics with genuine Kt signatures. For dijet events, the ocj variable is defined 
as: 



where Ej'^ is the transverse energy of the less energetic jet, and Mj is the transverse mass of 
the dijet system. For a perfectly measured dijet event with Ej^^ = ^t^^ arid jets back-to-back in 
(p, and in the limit in which each jet's momentum is large compared with its mass, the value of 
ax is 0.5. In the case of an imbalance in the measured transverse energies of back-to-back jets, 
oij is smaller than 0.5. Values significantly greater than 0.5 are observed when the two jets are 
not back-to-back, recoiling against genuine Kt- 

For events with three or more jets, an equivalent dijet system is formed by combining the jets 
in the event into two pseudo-jets. The Ej of each of the two pseudo-jets is calculated as the 
scalar sum of the measured Ej of the contributing jets. The combination chosen is the one 
that minimises the Ej difference (AHj) between the two pseudo-jets. This simple clustering 
criterion provides the best separation between multijet events and events with genuine Bj- 
Thus, in the case of events with at least three jets, the aj variable can be defined as: 

_ 1 Ht - AHt _ 1 1 - (AHt/Ht) 

Events with extremely rare but large stochastic fluctuations in the calorimetric measurements 
of jet energies can lead to values of ocj slightly above 0.5. Such events are rejected by requiring 
ax > 0.55. A similar behaviour is observed in events with reconstruction failures, severe energy 
losses due to detector inefficiencies, or jets below the £x threshold that result in significant fij 
relative to the value of ^ x (as measured by the calorimeter systems, which is not affected by 
jet £x thresholds). These classes of events are rejected by applying dedicated vetoes, described 
further in Ref. iT4| . The leakage above 0.5 becomes smaller with increasing Hj due to the 
increase in average jet energy and thus an improvement in jet energy resolution. Further, the 
relative impact of jets falling below the Ej threshold is reduced as the energy scale of the event 
(i.e. Hj) increases. 

The signal region is defined by Hj > 275 GeV and ax > 0.55, which is divided into eight bins 
in Hj: two bins of width 50 GeV in the range 275 < Hj < 375 GeV, five bins of width 100 GeV 
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5 Background estimation from data 



in the range 375 < Hj < 875 GeV, and a final open bin, Hj > 875 GeV. As in Ref. Cll, the jet 
Ej threshold is scaled for the two lowest Hj bins leading to thresholds of 37 GeV and 43 GeV. 
The two highest- Et jet thresholds are scaled to 73 GeV and 87 GeV. This approach maintains 
SM background admixtures and event kinematics similar to those observed for the higher Hj 
bins. Events are further categorised according to whether they contain exactly zero, one, two, 
or at least three reconstructed b-quark jets. 

Events in the signal region are recorded with a dedicated trigger condition that must satisfy si- 
multaneously the requirements Hj > 250 GeV and ccj > 0.53, with the latter threshold increas- 
ing to 0.60 towards the end of 2011 due to higher instantaneous luminosities. The efficiency 
with which events that would satisfy the signal region selection criteria also satisfy the trigger 
conditions is measured in data to be (82.8 ± 1.1)%, (95.9 ± 0.9)%, and (> 98.5 ± 0.9)% for the 
regions 275 < Hj < 325 GeV, 325 < Hj < 375 GeV, and Hj > 375 GeV, respectively 

A disjoint hadronic control sample consisting predominantly of multijet events is defined by 
inverting the aj requirement for a given Hj region, which is used primarily in the estimation 
of any residual background from multijet events. These events are recorded by a set of triggers 
with thresholds only in Hj. 

5 Background estimation from data 

Once all the signal region selection requirements have been imposed, the contribution from 
multijet events is expected to be negligible. The remaining significant backgrounds in the sig- 
nal region stem from SM processes with genuine in the final state. In the case of events 
where no b-quark jets are identified, the largest backgrounds with genuine U t arise from the 
production of W and Z bosons in association with jets. The weak decay Z — )■ is the only 
significant contribution from Z + jets events. For W + jets events, the two dominant sources are 
leptonic W decays in which the lepton is not reconstructed or fails the isolation or acceptance 
requirements, and the weak decay W — ?► rv where the t decays hadronically and is identified 
as a jet. Contributions from SM processes such as single-top, Drell-Yan, and diboson produc- 
tion are also expected. For events with one or more reconstructed b-quark jets, tt production 
followed by semi-leptonic weak decays becomes the most important single background source. 
For events with only one reconstructed b-quark jet, the contribution of both W + jets and Z + 
jets backgrounds are of a similar size to the tt background. For events with two reconstructed 
b-quark jets, tt production dominates, while events with three or more reconstructed b-quark 
jets originate almost exclusively from tt events, in which at least one jet is misidentified as 
originating from a bottom quark. 

In order to estimate the contributions from each of these backgrounds, three data control sam- 
ples are used, which are binned in the same way as the signal region. The irreducible back- 
ground of Z — > i/v + jets events in the signal region is estimated from two independent data 
samples of Z — ?> fij-i + jets and 7 + jets events, both of which share the kinematic properties 
ofZ — > 1/1/ + jets but have different acceptances. The Z ^ jifi + jets events have identical kine- 
matic properties to the Z — > vv + jets background when the two muons are ignored, but a 
smaller branching fraction, while the 7 + jets events have similar kinematic properties when 
the photon is ignored 1131 [381/ but a larger production cross section. A}i + jets data sample pro- 
vides an estimate for all other SM backgrounds, which is dominated by tt and W production 
leading to W + jets final states. 

The event selection criteria for the control samples are defined to ensure that any potential con- 
tamination from multijet events is negligible. Further, the same selection criteria also strongly 
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suppress contributions from a wide variety of SUSY models, including those considered in this 
analysis. Any potential signal contamination in the data control samples is accounted for in the 
fitting procedure described in Section |6] 

5.1 Definition of data control samples 

The fi + jets sample is recorded using two different trigger strategies, to account for evolving 
trigger conditions during the 2011 run. The hadronic trigger condition, combining Hj and aj, 
is used for the region 275 < Hj < 375 GeV. Here, the event selection, following closely the 
prescription described in Ref. | |39) , requires exactly one isolated muon that satisfies stringent 
quality criteria, with pj > 10 GeV and \f]\ < 2.1. In order for the trigger to be maximally 
efficient, the requirement ocj > 0.55 is also imposed. For the region Hj > 375 GeV, the trigger 
condition requires both a muon above a pj threshold as high as 40 GeV and Hj > 300 GeV. The 
muon must satisfy pj > 45 GeV in order for the trigger to be maximally efficient, at (91.3 ± 
0.1)%. The requirement ocj > 0.55 is again imposed when zero b-quark jets are reconstructed 
per event. For events in which at least one b-quark jet is reconstructed, no ocj requirement 
is used. This approach increases the statistical precision of predictions derived from event 
samples containing b-quark jets, while the impact of relaxing the aj requirement is tested with 
a dedicated set of closure tests described in Section |5!2l 

In addition to the requirements described above, further selection criteria are applied. The 
transverse mass of the muon and ^ j system must be larger than 30 GeV to ensure a sample rich 
in W bosons. The muon is required to be separated from the closest jet in the event by At] and 
A(p such that the distance AR = {Arj)^ + {A(p)^ > 0.5. To ensure that this sample is disjoint 
from the fifi + jets sample, the event is rejected if a second muon candidate is identified that 
does not satisfy all quality criteria or is non-isolated or is outside the acceptance, and the two 
muon candidates have an invariant mass that is within a window of ±25 GeV around the mass 
of the Z boson. 

The + jets sample follows the same trigger strategy and muon identification criteria as the 
}i + jets sample. The event selection requires exactly two oppositely charged, isolated muons 
satisfying stringent quality criteria, and an invariant mass within a window of ±25 GeV around 
the mass of the Z boson. Each muon is required to be separated from the nearest jet in the event 
by the distance AR > 0.5. The same aj requirements are used as for the fi + jets sample. 

The 7 + jets sample is selected using a dedicated photon trigger condition requiring a localised 
energy deposit in the ECAL with Ej > 135 GeV that satisfies loose photon identification and 
isolation criteria 1341 . The event selection requires Hj > 375 GeV, uj > 0.55, and a single 
photon to be reconstructed with Ex > 150 GeV, |//| < 1.45, satisfying tight isolation criteria, 
and with a minimum distance to any jet of AR > 1.0. For these selection criteria, the photon 
trigger condition is found to be fully efficient. 

5.2 Method and systematic studies 

The method used to estimate the SM background contributions in the signal region relies on the 
use of translation factors, which are functions of Hj and the number of b-quark jets per event, 
Mb, and are computed separately for each data control sample. These translation factors are de- 
termined from simulation samples generated with MadGraph v4.22 1401 interfaced to pythia 
6.4 tune Z2 IHI, and the GEANT 4-based ||43 CMS detector simulation. Each factor is defined 
as the ratio of yields from simulation in a given bin of the signal region, N^^^^{Hj, n\,) and the 
corresponding bin of one control sample, N^Q^°^{Hj,n^,). The factors are used to translate the 
observed yield measured in a control sample bin, N^^^^^°^{Hj, n^) into an expectation for one or 
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5 Background estimation from data 



more SM background processes in the corresponding bin of the signal region, N^^^^{Hj, n\,): 

signal. „ N _ .rcontrol/TT „ \ ^ ^MC i^T'^h) 

In order to maximise sensitivity to potential new physics signatures in final states with multiple 
b-quark jets, a method that improves the statistical power of the predictions from simulation, 
particularly for > 2, is employed. The distribution of is estimated from generator-level in- 
formation contained in the simulation, namely the number of reconstruction-level jets matched 
to underlying b quarks, n^^'^, and light quarks, n^^^, per event. All relevant combinations of 
n^^'^ and n^'^ are considered, and event counts are recorded in bins of Hj for each combination 
N [n^^^ , fiq'^ , Hj) . The b-tagging efficiency, e, and a flavour-averaged mistagging rate, m, are 
also determined from simulation for each Hj bin, with both quantities averaged over jet pj 
and f]. Corrections are applied on a jet-by-jet basis to both e and m in order to match the cor- 
responding measurements with data 15511361 . This information is sufficient to predict and 
thus also determine the yield from simulation for a given bin, N{Hj, n^,): 



N(Ht, nb) = E N(nf ^ Hj) x P{nl'S. „gen^ ^ p^^tag. ^gen^ 

where n^^ and n^q^ are the number of times a reconstruction-level b-quark jet originates from 
an underlying b-quark and light-quark respectively, and P(n^^^; n^^", e) and P(nq ^; n^", m) are 
the binomial probabilities for this to happen. The predicted yields are found to be in good 
agreement with the yields obtained directly from the simulation in those bins with significant 
population. The method exploits the ability to make precise measurements of N{n^^, nq*^'^, Hj), 
e, and m independently of n^, which means that event yields for a given b-quark jet multiplicity 
can be predicted with a higher statistical precision than obtained directly from simulation. A 
precise measurement of m is particularly important for events with > 3, which occurs in the 
SM because of the presence of mistagged jets in the event. In this case, the largest background 
is tt, with two correctly tagged b-quark jets and an additional mistagged jet. 

Since the translation factors are obtained from simulation, an appropriate systematic uncer- 
tainty is assigned to each factor to account for theoretical uncertainties [38J and limitations in 
the simulation modelling of event kinematics [13J. The magnitudes of the systematic uncer- 
tainties are determined from a representative set of closure tests in data, in which yields from 
one of the three independent control samples, along with the corresponding translation factors 
obtained from simulation, are used to predict the yields in another control sample, following 
the same prescription defined in Eq. (j3|. 

A set of five closure tests use the three data control samples to probe key ingredients of the 
simulation modelling of the SM backgrounds with genuine £ t as a function of Hj, as shown 
in Fig. [T] The first three closure tests are carried out within the }i + jets sample, and probe 
the modelling of the ocj distribution in genuine £ j events (circles), the relative contributions 
of W + jets and tt events (squares), and the modelling of the reconstruction of b-quark jets 
(triangles), respectively. The fourth test (crosses), connecting the }i + jets and }i}i + jets control 
samples, addresses the modelling of the relative contributions of Z + jets to the sum of both W 
+ jets and tt events, while the fifth test (stars) deals with the consistency between the Z— > 
+ jets and 7 + jets samples. All individual closure tests demonstrate, within the statistical 



7 



T3 
Q. 



to 
o 



CMS, L =4.98fb"^ 



-I— p-i— I— I— I— p 
\^ = 7 TeV 



Q Systematic uncertainty 
O < 0.55 -> otj > 0.55 (|i + jets) 
n b tags 1 b tag + jets) 
A 1 b tag 2 b tags (i^ + jets) 
O |j. + jets -> |in + jets 
■jir HH + jets ^ y + jets 



.f^ i^n ,1^11 



_L 



_L 



300 400 500 600 700 800 900 

Hj (GeV) 

Figure 1: A set of five closure tests, described in the text, that use the three data control samples 
to probe key ingredients of the simulation modelling of the SM backgrounds, as a function 
of Hj. Error bars represent statistical uncertainties only. The shaded bands represent the Hj- 
dependent systematic uncertainties assigned to the translation factors. 



precision of each test, that there are no significant biases inherent in the translation factors 
obtained from simulation. The level of closure achieved in these tests is used to estimate the 

systematic uncertainties that are assigned to the translation factors, which are determined for 
three regions 275 < Hj < 575 GeV, 575 < Hj < 775 GeV, and Hj > 775 GeV to be 10%, 20%, 
and 40%, respectively. 

A further dedicated study to account for potential systematic effects arising from the modelling 
of the reconstruction of b-quark jets in the simulation has been performed. After correcting the 
efficiency and mistagging rates of b-quark jets in simulation for residual differences as mea- 
sured in data, the corresponding uncertainties on these corrections are propagated to the trans- 
lation factors and found to be at the sub-percent level. In addition, several robustness tests are 
performed, including treating c-quark jets as b-quark jets in the yield estimates throughout, 
as well as ignoring the contribution from hadronic r-lepton decays. These tests also demon- 
strate sub-percent effects on the translation factors, highlighting the insensitivity to potential 
mismodelling in simulation. Hence, the Hj-dependent systematic uncertainties of 10%, 20%, 
and 40% are used for all b-quark jet multiplicities. 



6 Results 

A likelihood model of the observations in all four data samples is used to obtain a consistent 
prediction of the SM background, and to test for the presence of a variety of signal models. It 
is written as 

Ltotal = rr ( i-hadronic ^ ^ Ajets ^ ^m+iets ^ ^7+jets ) 

«b=0 



^+jets ii}i+)ets 7+jets J 
hadronic ^ ^^+jets 



^ ^hadronic ^ ^w+iets ' (^) 
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6 Results 



where i-hadronic describes the yields in the eight Hj bins of the signal region when exactly n\, 
reconstructed b-quark jets are required. In each bin of Hj, the observation is modelled as 
Poisson-distributed about the sum of a SM expectation and a potential signal contribution. The 
components of this SM expectation are related to the expected yields in the control samples via 



translation factors derived from simulation, as described in Section 5.2 Signal contributions in 
each of the four data samples are considered, though the only significant contribution occurs 
in the signal region and not the control samples. The systematic uncertainties associated with 
the translations are accounted for with nuisance parameters, the measurements of which are 
treated as normally-distributed. Since for tib > 3 the dominant SM background arises from 
top events, only the }i + jets control sample is used in the likelihood to determine the total 
contribution from all (non-multijet) SM backgrounds in the signal region. 

In addition, any potential contribution from multijet background in the signal region is ac- 
counted for by using the ratio of events which result in a value of ocj above and below some 
threshold value for a given Hj bin. The dependence of this ratio, R^j, on Hj is modelled as a 
falling exponential function: A^^e^'^^T |[14|. A common parameter k is used for all four cate- 
gories of b-quark jet multiplicity, and is constrained via measurements in a multijet-enriched 
data side-band satisfying the criteria Hj < 575 GeV and 0.52 < uj < 0.55. A further side-band, 
defined by inverting the Ht/^t requirement of Ref. [14J, is used to confirm that this method 
provides an unbiased estimate of k and to determine a systematic uncertainty. 

In order to test the compatibility of the observed yields with the expectations from SM pro- 
cesses only, signal contributions are fixed to zero and the likelihood function is maximised 
over all parameters. The maximum likelihood values of the multijet normalisation parameters 
are found to be compatible with zero, within uncertainties, confirming the hypothesis that 
the multijet background is negligible after the final selection. Further, the SM expected yields 
obtained from an alternate fit, in which these normalisation parameters are fixed to zero, agree 
well with those obtained from the nominal fit. 

The signal region data yields, as well as the SM expectations obtained from the simultaneous 
fit across all samples, are shown in Table [l] A comparison of the observed yields and the SM 
expectations in bins of Hj for events with exactly zero, one, two, and at least three reconstructed 
b-quark jets are shown in Figures |2j|3||4| and|5j respectively, for the signal region and the three 
control samples. In all four categories of b-quark jet multiplicity, the samples are well described 
by the SM hypothesis. In particular, no significant excess above the SM expectation is observed 
in the signal region. 

Table 1: Comparison of the observed yields in the different Hj and b-quark jet multiplicity 
bins for the signal region with the SM expectations and combined statistical and systematic 
uncertainties given by the simultaneous fit. 
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Figure 2: Comparison of the observed yields and SM expectations given by the simultaneous 
fit in bins of Hj for the (a) signal region, (b) fi + jets, (c) fiji + jets, and (d) 7 + jets samples when 
requiring exactly zero reconstructed b-quark jets. The observed event yields in data (black dots) 
and the expectations and their uncertainties, as determined by the simultaneous fit, for all SM 
processes (light blue solid line with dark blue bands) are shown. For illustrative purposes only, 
the signal expectation (magenta dashed line) in the signal region for the simplified model A 
(defined in Section [7!2j | with nig = 800 GeV and wJLgp = 200 GeV is superimposed on the SM 
background expectation. 
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Figure 3: Comparison of the observed yields and SM expectations given by the simultaneous 
fit in bins of Hj for the (a) signal region, (b) + jets, (c) + jets, and (d) 7 + jets samples. 
Same as Fig. |2| except requiring exactly one reconstructed b-quark jet. The observed event 
yields in data (black dots) and the expectations and their uncertainties, as determined by the 
simultaneous fit, for all SM processes (light blue solid line with dark blue bands) are shown. 
For illustrative purposes only, the signal expectation (magenta solid line) in the signal region 
for the simplified model D (defined in Section [7!2] | with nig = 500 GeV and Wlsp = 150 GeV is 
superimposed on the SM background expectation. 
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Figure 4: Comparison of the observed yields and SM expectations given by the simultaneous 
fit in bins of Hj for the (a) signal region, (b) }i + jets, (c) + jets, and (d) 7 + jets samples. 
Same as Fig. |2j except requiring exactly two reconstructed b-quark jets. The observed event 
yields in data (black dots) and the expectations and their uncertainties, as determined by the 
simultaneous fit, for all SM processes (light blue solid line with dark blue bands) are shown. 
For illustrative purposes only, the signal expectation (magenta solid line) in the signal region 
for the simplified model D (defined in Section [7!2| | with nig = 500 GeV and Wlsp = 150 GeV is 
superimposed on the SM background expectation. 
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Figure 5: Comparison of the observed yields and SM expectations given by the simultaneous 
fit in bins of Hj for the (a) signal region and (b) fi + jets samples. Same as Fig.|2j except requiring 
at least three reconstructed b-quark jets. The observed event yields in data (black dots) and the 
expectations and their uncertainties, as determined by the simultaneous fit, for all SM processes 
(light blue solid line with dark blue bands) are shown. For illustrative purposes only, the signal 
expectation (magenta solid line) in the signal region for the simplified model D (defined in 
Section 7.2 1 with nig = 500 GeV and jwlsp = 150 GeV is superimposed on the SM background 
expectation. 



7 Interpretation of the results 

Limits are set in the parameter space of the CMSSM and in a set of simplified models that char- 
acterise both third-generation squark production and compressed SUSY spectra scenarios. The 
CLs method [43, 44J is used to compute the limits, with the one-sided profile likelihood ratio 
as the test statistic 1451 . The sampling distributions for the test statistic are built by generating 
pseudo-data from the likelihood function, using the respective maximum-likelihood values of 
the nuisance parameters under the background-only and signal-plus-background hypotheses. 

Events samples for the CMSSM and simplified models are generated at leading order with 
PYTHIA 6.4 [|4T|. Inclusive, process-dependent, next-to-leading order calculations with next- 
to-leading logarithmic corrections Il46ti50l (NLO+NLL) of SUSY production cross sections are 
obtained with the program PROSPINO fST] and CTEQ6M [52J parton distribution functions. The 
simulated signal events include multiple interactions per LHC bunch crossing (pileup) with the 
distribution of reconstructed vertices that match the one observed in data. 

7.1 Interpretation in the CMSSM 

The CMSSM is described by the following five parameters: the universal scalar and gaugino 
mass parameters, niQ and mi/2; the universal trilinear soft SUSY-breaking parameter, Aq; the ra- 
tio of the vacuum expectation values of the two Higgs doublets, tan /6; and the sign of the Higgs 
mixing parameter, }i. At each point in the parameter space of the CMSSM, the SUSY particle 
spectrum is calculated with SoftSUSY I|53| . Experimental uncertainties on the SM background 
prediction (10-40%), the luminosity measurement (2.2%) [54|, and the total selection efficiency 
times acceptance for the considered signal model (16%) are included in the calculation of the 
limit. The dominant sources of uncertainty on the signal efficiency times acceptance are de- 
rived from systematic variations of parton distribution functions, and corrections applied to jet 
energies and b-tagging efficiency and mistag rates. 



7.2 Interpretation with simplified models 
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Figure |6] shows the observed and expected exclusion limits at 95% confidence level (CL) in the 
(mo,mi/2) plane for tan ,6 = 10 and Aq = OGeV, calculated with the NLO+NLL SUSY pro- 
duction cross section. For this choice of parameter values, squark masses below 1250 GeV are 
excluded at 95% CL, as are gluino masses below the same value for the region mo < 600 GeV. In 
the region 600 < niQ < 3000 GeV, gluino masses below 700 GeV are excluded, while the squark 
mass in the excluded models varies in the range 1250-2500 GeV, depending on the value of 
niQ. The mass limits are determined conservatively from the observed exclusion based on the 
theoretical production cross section minus \cr uncertainty [|55|. 
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Figure 6: Exclusion contours at 95% CL in the CMSSM (mo, mi/2) plane (tan/3 = 10, Aq = 
Qfji > Q) calculated with NLO+NLL SUSY production cross sections and the CLg method. The 
solid black line indicates the observed exclusion region. The dotted-dashed black lines repre- 
sent the observed excluded region when varying the production cross section by its theoretical 
uncertainty. The expected median exclusion region (green dashed line) ±lc7" (green band) are 
also shown. The CMSSM template is taken from Ref. Ii56i . 



7.2 Interpretation with simplified models 

The data observations are also interpreted using simplified models that characterise third- 
generation squark production and compressed spectra scenarios, where the mass difference 
between the primary produced sparticle (e.g. a squark or a gluino) and the LSP is rather small. 
The production and decay modes of the models under consideration are summarised in Ta- 
ble |2] The simplified models A and B are used to characterise the pair production of gluinos 
and first- or second-generation squarks, respectively, depending on their mass as well as on the 
LSP mass. Simplified models C to f describe various production and decay mechanisms in the 
context of third-generation squarks. 

Experimental uncertainties on the SM background predictions (10^0)%, the luminosity mea- 
surement (2.2%), and the total acceptance times efficiency of the selection for the considered 
signal model (12%— 18%) are included in the calculation of the limit. Signal efficiency in the 
kinematic region defined by < »Zg(q) — Wlsp < 175 GeV or Mtg(q) < 300 GeV is due in part 
to the presence of initial-state radiation. Given the large associated uncertainties, no interpre- 
tation is provided for this kinematic region. In the case of model E, for which pair-produced 
gluinos decay to tt pairs and the LSP, the region < mg — niigp < 400 GeV is not considered. 
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Figure 7: Upper limit on cross section at 95% CL as a function of mq or mg and OTlsp for various 
simplified models. The solid thick black line indicates the observed exclusion region assum- 
ing NLO+NLL SUSY production cross section. The thin black lines represent the observed 
excluded region when varying the cross section by its theoretical uncertainty. The dashed pur- 
ple lines indicate the median (thick line) ±lcr (thin lines) expected exclusion regions. The mass 
ranges considered for models C and E differ from the other models. 
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Table 2: The first three columns define the production and decay modes for various simplified 
models. The last two columns indicate the search sensitivity for these models, where 

and represent the largest mass beyond which no limit can be set for squarks/ gluinos and 
the LSP, respectively. The exclusion range for wiqj-g) is bounded from below by the kinematic 
region considered for each simplified model, as defined in the text. The quoted estimates are 
determined conservatively from the observed exclusion based on the theoretical production 
cross section minus Icr uncertainty. For model C, the search is at the threshold of sensitivity for 
the considered (mq, wilsp) parameter space, as discussed in the text. 
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Figure [7| shows the upper limit on the cross section at 95% CL as a function of mq or mg 



and 



Wlsp for various simplified models. The point-to-point fluctuations are due to the finite number 
of pseudo-experiments used to determine the observed upper limit. The solid thick black line 
indicates the observed exclusion region assuming NLO+NLL SUSY cross section for squark 
pair production in the limit of very massive gluinos (or vice versa). The thin black lines repre- 
sent the observed excluded region when varying the cross section by its theoretical uncertainty. 
The dashed purple lines indicate the median (thick line) ±1 c (thin lines) expected exclusion 
regions. 

The most stringent mass limits on the pair-produced sparticles are obtained at low LSP masses, 
while the limits typically weaken for compressed spectra, i.e., points close to the diagonal. In 
particular, for all of the considered simplified models, there is an LSP mass beyond which no 



limit can be set. This is illustrated in Figure 7a where the most stringent limit on the gluino 
mass is obtained at around 950 GeV for low LSP masses, while this limit weakens to below 
900 GeV when the LSP mass reaches 350 GeV. For LSP masses above 400 GeV, no gluino masses 
can be excluded. Tabled summarises these two extreme cases for models A to F. The estimates 
on the mass limits are determined conservatively from the observed exclusion based on the 
theoretical production cross section minus la uncertainty. 

No exclusion of direct top squark pair production (model C) assuming the NLO+NLL pro- 
duction cross section is expected with the analysed dataset and for LSP masses greater than 
50 GeV. Figure |8] shows the observed upper limit at 95% CL on the cross section as a function 
of the top squark mass (m^) only, for a fixed LSP mass of nii^sp = 50 GeV. Within the mass range 
350 < wjj < 475 GeV, the observed upper limit fluctuates about the theoretical production cross 
section minus la uncertainty. This mass range is fully excluded when considering the nominal 
production cross section. 



8 Summary 

A search for supersymmetry using the CMS detector is reported, based on a data sample of 
PP collisions collected at ^/s = 7TeV, corresponding to an integrated luminosity of 4.98 ± 
0.11 fb^^. Final states with two or more jets and significant Bt, as expected from high-mass 
squark and gluino production and decays, have been analysed. An exclusive search has been 
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Figure 8: Excluded cross section versus top squark mass for a model in which pair-produced 
top squarks decay to two top quarks and two neutralinos of mass wilsp = 50 GeV. The solid 
blue line indicates the observed cross section upper limit (95% CL) as a function of the top 
squark mass, m^. The dashed orange line and blue band indicate the median expected excluded 
cross section with experimental uncertainties. The solid black line and grey band indicate the 
NLO+NLL SUSY top squark pair-production cross section and theoretical uncertainties. 

performed in a binned signal region defined by the scalar sum of the transverse energy of jets, 
Hj, and the number of jets identified to originate from a bottom quark. The sum of standard 
model backgrounds per bin has been estimated from a simultaneous binned likelihood fit to 
hadronic, }i + jets, ]A}i + jets, and 7 + jets samples. The observed yields are found to be in agree- 
ment with the expected contributions from standard model processes. Limits in the CMSSM 
(mo, mi/2) plane for tanj6 = 10, Aq = OGeV, and ji > have been derived. For this choice 
of parameter values, gluino masses below 700 GeV are excluded at 95% CL. The exclusion in- 
creases to 1250 GeV for squarks and gluinos of comparable mass. Furthermore, exclusion limits 
are also set in simplified models, with a special emphasis on third generation squarks and com- 
pressed spectra scenarios. In the considered models with gluino pair production and for small 
LSP masses, typical exclusion limits of the gluino mass are around 1 TeV. For simplified mod- 
els with squark pair production, first or second generation squarks are excluded up to around 
750 GeV and bottom squarks are excluded up to around 500 GeV, again for small LSP masses. 
No exclusion is expected for direct pair production of top squarks that each decay to a top quark 
and a neutralino of mass niisp > 50 GeV. However, within the mass range 350 < m ^ < 475 GeV 
and for Wlsp = 50 GeV, the observed upper limit fluctuates about the theoretical production 
cross section minus 1 a uncertainty. Thus, for the simplified models under consideration, the 
most constraining limits on the LSP and third-generation squark masses indicate that a large 
range of SUSY parameter space is yet to be probed by the LHC. 
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